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It is generally agreed that the bursting phenomenon is the most significant 
dynamic event in a turbulent boundary layer. About 80 percent of the momentum 
transport occurs during these bursts. Previous attempts to understand the 
physics and structure of these events were frustrated by the fact that bursts 
occur randomly in space and time and that successive bursts are not 
necessarily identical. During Phase I of this investigation, "artificial" 
bursts were generated in laminar and turbulent boundary layers. The 
burst-lilce events were produced by withdrawing near-wall fluid from two minute 
holes separated in the spanwise direction or by pitching a miniature delta 
wing that was flush-mounted to the wall. Either of these actions generated 
streamwise vorticity and a low-speed streak that resembled a naturally 
occurring one. The resulting sequence of events occurred at a given location 
and at controlled times, allowing detailed examination and comparison with 
natural, random bursts by means of flow visualization and fast-response probe 
measurement techniques. 

By all available measures, the Phase I research should be considered highly 
successful. We have generated an artificial burst in a turbulent boundary 
layer. This will allow detailed study of the dynamics of a bursting event and 
make it easier to develop methods to control or eliminate bursts, thus 
reducing drag. We are currently contemplating a novel method to achieve 
substantial drag reduction. The device combines the beneficial effects of 
suction and a longitudinally ribbed surface, and will be optimized using 
artificially-generated bursts. A successful application of the proposed 
innovation on commercial aircraft will result in an annual fuel savings of 
several billion dollars. 
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1 . Introduction 

Recent turbulent boundary layer research has clearly shown that the wall 
region is dominated by a sequence of eddy motions that are collectively called 
the bursting phenomenon. This process was reviewed by Willmarth (1975) and 
summarized by Blackwelder (1978). Qualitatively, the process begins with a 
pair of elongated, streamwise, counter-rotating vortices having diameters of 
approximately 40 v/u^, where V/u^ is the viscous scale, v is the kinematic 
viscosity and u^ is the friction velocity. The vortices exist in a strong 
shear and induce low- and high-speed regions between them as shown in 
Section BB of Figure 1. The vortices and the accompanying eddy structures 
occur randomly in space and time. However, their appearance is regular 
enough that an average spanwise wavelength of approximately 80 to 100 v/u^. has 
been identified by Kline et al. (1967) and others. Kline et al. also observed 
that the low-speed regions grow downstream and develop inflectional U(y) pro- 
files, as sketched in Figure 1. At approximately the same time, the interface 
between the low- and high-speed fluid begins to oscillate. The low-speed 
region lifts up away from the wall as the oscillation amplitude increases, and 
then the flow rapidly breaks down into a completely random pattern. Since 
this latter process occurs on a very short time scale, Kline et al. called it 
a "burst." Corino & Brodkey (1969) showed that the low-speed regions are 
quite narrow, i.e., 20 v/u , and may also have significant shear in the 
spanwise direction. 

Considerably more has been learned about the bursting process during the 
last decade. For example, Falco (1980) has shown that the wall region is con- 
tinuously bombarded by "pockets" of high-speed fluid originating in the loga- 
rithmic and possibly the outer layers of the flow. These pockets tend to 
promote and/or enhance the inflectional velocity profiles by increasing the 
instantaneous shear leading to a more rapidly growing instability. Blackwelder 
& Haritonidis (1983) have shown that the frequency of occurrence of these 
events scales with the viscous parameters consistent with the usual boundary 
layer scaling arguments. 

The complexity of the bursting phenomenon should be apparent from the fact 
that extensive investigations have failed to provide an explanation for the 
origin of bursting (Offen & Kline, 1975; Praturi & Brodkey, 1978; Falco, 1979; 
1980; Smith, 1983) or the topological characteristics of the bursting coherent 
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Figure 1. Sketch of the Low-Speed Streaks in the Wall Region. The 

Counter Rotating Eddies between the Streaks are Indicated in 
Section B-B. 
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structure, as well as its evolution and dynamics (Hussain, 1983a). Sophisti- 
cated experiments have been conducted to capture the bursting signature 
(Blackwelder & Kaplan, 1976; Nishioka et al., 1981; Antonia, 1981), but it is 
clear that the educed signature is dependent on the detection scheme as well 
as on the threshold level of a given detection scheme (Kunen et al., 1983; 

Sato, 1983). 

Hussain (1983b) reiterated that a single-point detection of an inherently 
three-dimensional random structure must produce significant smearing. That 
is, eduction of the bursting coherent structure should not only involve a 
three-dimensional detection scheme but also require first the identification 
of the mode (i.e., spatial configuration) of the structure and then the pro- 
bability density distribution as a function of the characteristic structure 
parameters, like shape, size, strength, orientation, convection velocity, etc. 
Furthermore, a stationary sensor may intersect successive three-dimensional 
structures (say, hairpin vortices) at different relative spanwise locations, 
and a structure may advect past a sensor at different transverse displacements 
so that a stronger structure farther away could not be necessarily differen- 
tiated from a weaker structure nearby. 

It should thus be clear that an array of probes in a transverse rake must 
be used for measurement. Furthermore, the measuring rake should be located in 
between two detection rakes so that the relative spanwise location of the 
measuring rake within the three-dimensional advecting structure can be deter- 
mined instantaneously. 

Since the coherent vorticity 0 is the principal property of a coherent 
structure, structure detection must be based on and not on the correlation 
of velocity, pressure and/or intermittency signals (Hussain 1980). Only large- 
scale vorticity is involved, however, so a transverse rake of X-wires is 
adequate to measure the instantaneous map of the spanwise vorticity <fl> in the 
(y, t) plane. From the smoothed maps <fi(y, t)>, one should decide which struc- 
tures should be accepted. That is, eduction involves accepting only those 
structures of a given mode with the peak value (denoting strength) of <ft> above 
a certain level and the peak location at an optimally selected transverse loca- 
tion and accepting structures of a given shape and size. The accepted struc- 
tures should then be ensemble averaged after proper phase alignment. One way 
to align is appropriate time shifting of realizations via the method of 
iterative cross-correlation of each realization with the ensemble average 
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(Zilberman et al., 1977). Whatever survives the phase averaging (i.e. , 
ensemble averaging at a fixed phase) of the appropriately aligned realizations 
is the coherent structure, and the departure of each realization from the 
phase averaging is the corresponding incoherent turbulence . The phase- 
averaged contours of coherent vorticity, coherent and incoherent Reynolds 
stress, coherent production, etc., provide the necessary structure properties. 
The phase-averaged velocity vector pattern in the frame of the structure super- 
imposed on the contours of <£2> will denote the spatial variations of entrain- 
ment (Hussain & Zaman, 1980). 

It is necessary to emphasize that these elaborate efforts are justified if 
the flow is characterized by "preferred modes", as the coherent structure 
approach to turbulence is not very helpful if the flow has a large variety of 
coherent structures. Furthermore, capturing an individual structure via 
spatial smoothing is not helpful, as there is no test for convergence of the 
smoothing that removes sharp fronts otherwise retained by ensemble averaging. 
Also, one must assure that the educed structure is not "freak" and does indeed 
represent a "preferred mode". A classical example of such a freak structure 
is the single, large horseshoe vortex detected by Wygnanski et al. (1976) and 
Cantwell et al. (1978) by ensemble averaging an artificial turbulent spot. The 
results of Gad-el-Hak et al. (1981) clearly indicate the dynamic insignificance 
of such a structure. 

Returning to the bursting coherent structure, it should be clear that the 
variations in their characteristic parameters as well as the random occurrence 
in time and the random transverse and spanwise displacements of natural 
bursting structures must have produced unacceptable smearing, thus frustrating 
prior eduction schemes in the boundary layer. Furthermore, we are unaware of 
any attempt to determine if all structures included in the educed bursting 
signature (which differed from a phase average) are of the same mode. In an 
attempt to remove significantly the inherent smearing in the eduction of 
natural bursting structures, a new approach has been attempted: eduction of 

artificially induced bursting . The structure details can then be educed via 
measurements phase-locked to the excitation (Hussain & Reynolds, 1972), much 
the same as was done in the case of the artificial turbulent spot (Gad-el-Hak 
et al., 1981). During the first phase of this work, we were able to generate 
an artificial burst in both laminar and turbulent boundary layers. The 
sequence of events induced artificially was compared to that of a natural 
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burst in a turbulent boundary layer using flow visualization and hot-film 
probe measurements. In all the comparisons that we have made thus far, the 
artificial event was kinematically and dynamically identical to the bursts 
occurring randomly in space and time in a turbulent boundary layer. 
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2 . Experimental Approach 

2.1 Towing Tank 

Turbulent and laminar boundary layers were generated by towing a flat plate 
in a water channel that is 18 meters long, 1.2 meters wide and 0.9 meter deep. 
The towing tank, shown in Figure 2, has been described by Gad-el-Hak et al. 
(1981). The flat plate was rigidly mounted under a carriage that rides on two 
tracks mounted on top of the tank. During towing, the carriage was supported 
by an oil film to ensure a vibrationless tow, having an equivalent freestream 
turbulence of about 0.1 percent. The carriage was towed by two cables driven 
through a reduction gear by a 1.5-hp Boston Ratiotrol motor. The towing speed 
was regulated within an accuracy of 0.1 percent. The system was able to 
achieve towing speeds between 5 and 140 cm/sec for the present study. However, 
most of the runs reported here were conducted at a speed of 20 cm/sec. 

2. 2 Flat Plate 

A unique, zero-pressure-gradient flat plate was constructed for the present 
investigation. The 1-meter by 2-meter structure, shown in Figure 3, is made 
of glass-reinforced polycarbonate plate 6 mm thick, glued to a stainless steel 
frame designed for minimum obstruction to the flow. The plate has an elliptic 
nose at the leading edge and an adjustable lifting flap at the trailing edge. 

To avoid leading-edge separation and premature transition, the flap is adjusted 
so that the stagnation line near the leading edge is located on the working 
surface of the plate. The plate's smooth surface and flatness to within a few 
microns make it one of the most controlled test beds available for boundary 
layer research. A laminar boundary layer can be obtained over the entire 
working surface for towing speeds in the range of 5 to 80 cm/sec. Trips are 
used to generate a fully developed turbulent boundary layer. The trips are 
brass cylinders 0.32 cm in diameter and 0.25 cm high placed 5 cm downstream of 
the leading edge and having their axes perpendicular to the flat plate. 

2.3 Excitation Techniques 

To generate an artificial burst, several excitation techniques were tried. 
Suction was applied suddenly through one, two, three or four holes separated 
in the spanwise direction by a distance of AZ + = 100. The minute holes 
were 0.4 mm in diameter and were connected to a vacuum chamber controlled with 
a solenoid valve that is driven by a signal generator. This allowed the sudden 
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withdrawal of a given amount of fluid and the generation of a horseshoe vortex 
that evolved into a burst. 

The second excitation technique tried was the sudden pitching of a minia- 
ture delta wing having a span of 80 wall units. The wing, with its apex facing 
the flow, was normally flush-mounted to the surface of the plate. At the 
desired instant, a pin centered around the trailing edge of the wing was driven 
upward using a solenoid to pitch the delta wing to a negative angle of attack 
of 30°. This produced two strong, longitudinal, counter-rotating vortices that 
lifted away from the wall and evolved into a burst. The burst generators were 
located 79 cm downstream of the plate's leading edge. Figure 4 is a close-up 
photograph of the surface of the flat plate, showing the flush-mounted delta 
wing and the four suction holes (located to the left of the white marks in the 
picture) . 

2.4 Visualization Methods 

To visualize the bursting events in laminar and turbulent boundary layers, 
both fluorescent dye and hydrogen-bubble techniques were used. The dyes were 
seeped into the boundary layer through a spanwise slot 0.15 mm wide and 15 cm 
long and located 77 cm downstream of the leading edge of the flat plate. The 
slot is shown in the photograph in Figure 4 just upstream of the burst genera- 
tors. The dye slot was milled at a 45° angle inclined toward the plate's 
trailing edge to minimize the flow disturbance. The hydrogen bubbles were 
generated using a stainless steel wire having a diameter of 50 microns and a 
length of 10 cm. The wire was placed at different locations either parallel 
to or perpendicular to the flat plate to obtain top views or side views, re- 
spectively, of the flow field. To create time lines, a standard circuit was 
employed to supply intermittently 100 volts to the wire, which acted as a 
cathode. 

Both the fluorescent dyes and the hydrogen bubbles were illuminated using 
a sheet of laser light projected in the desired plane. This provided an extra 
degree of freedom in observing the bursting events, because both the tracer and 
the light location could be controlled within the limitations of the experi- 
mental apparatus. In order to generate a sheet of light, a 5-watt argon-ion 
laser (Spectra Physics Model 164) was used with a mirror mounted on an optical 
scanner having a natural frequency of 720 Hz (General Scanning Inc.) and driven 
by a sine-wave signal generator of the desired frequency. The frequency of the 
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Figure 4. Close Up of the Flush-Mounted Delta Wing 
and the Suction Holes 
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sine wave was usually set equal to the inverse of the shutter speed of the 
camera. The light sheets were approximately 1 mm thick, which was sufficient 
to resolve the large structure within the turbulent regions. A vertical sheet 
of laser light parallel to the flow was used to visualize side views of the 
bursting event, and a horizontal sheet near the surface of the plate was used 
to obtain top views of the flow field. 

2.5 Velocity Measurements 

Miniature boundary layer hot-film probes (TSI Model 1260) were used in the 
present investigation to measure the longitudinal mean and fluctuating velo- 
cities. The probe diameters were 0.025 mm, and their sensing lengths were 
0.25 mm. A probe traverse powered by a stepping motor and controlled through 
an APPLE II microcomputer was used for surveying the boundary layers. Conven- 
tional statistical quantities, such as the mean and the root mean square, were 
computed from the velocity signals. 

2.6 Pattern Recognition Algorithm s 

A new pattern recognition technique that utilizes the streamwise velocity 
signal from three hot-film probes was developed to detect low-speed streaks 
(Gad-el-Hak et al. , 1984). The probes were located at y + = 10 at the same 
streamwise position with a spanwise separation of approximately 20 v/u^. Thus, 
the total array spanned 40 V/u^, which is approximately half the average low- 
speed streak spacing. A pattern recognition algorithm identified a low-speed 
region whenever all three signals were less than the local mean velocity by 
one-half the root-mean-square (rms) value, and the velocity measured by the 
middle probe was less than the velocity on either side. The beginning and end 
of each low-speed region could thus be identified, and the streaks' average 
length and frequency computed. 

To test this algorithm, experiments were first carried out using a rake of 
twelve velocity probes located at y + = 10 but having a significantly broader 
spanwise extent. Low-speed streaks were readily identified using the isovelo- 
city contours computed from the rake's output, and their spacing was in reason- 
able agreement with results obtained from visualization experiments conducted 
under similar conditions. The results of the algorithm were compared with the 
low-speed streaks obtained from the rake. The comparison showed that the 
detected algorithm picked out the streaks exceedingly well, and that the 
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streaks meander in the spanwise direction so that often a single streak 
crosses the probes two or more times, resulting in multiple detections. 

The second pattern-recognition algorithm employed in the present investiga- 
tion was a burst detection scheme using the variable-interval time-averaging 
(VITA) technique developed by Blackwelder & Kaplan (1976). A single hot-film 
probe located at y + = 20 was used for burst detection. The program counted 
the number of bursts that occur near the wall and recorded their intensities. 
The low-speed streak algorithm and the burst detection program are reproduced 
in Appendix I. The FORTRAN listing also includes a conventional statistics 
program whose output is used by the burst and streak detectors. 
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3. Bursts in Laminar Flow 


The aim of the present investigation is to generate artificially a bursting 
event in a turbulent boundary layer. The resulting sequence of events occurs 
at a given location and at controlled times, thus allowing detailed examina- 
tion. To understand how the artificial burst generator works, the initial ex- 
periments were conducted in a laminar boundary layer environment. This allowed 
relatively easy visualization of the artificial events, free from the random 
background that exists in a turbulent boundary layer. Of course, the artifi- 
cial event evolves differently in the laminar and the turbulent cases. 

However, the resulting patterns are qualitatively the same. 

In a two-dimensional laminar boundary layer, the vorticity vector is in 
the spanwise direction. Consider the boundary layer to consist of a row of 
rectilinear vortices aligned in the spanwise direction. Consider further the 
result of suddenly withdrawing near-wall fluid from one, two or three suction 
holes separated in the spanwise direction. As shown in the schematic in 
Figure 5, a near-wall vortex line will then be "pinned down" at the location 
of the suction holes, thus forming one or more horseshoe-type hoops. In the 
case of two holes, one horseshoe vortex with its head pointing downstream 
forms. For three suction holes, two such vortices are generated. The longi- 
tudinal vortices exist in a strong shear and induce low- and high-speed regions 
between them, as shown in Figure 5. For the single-hole case, two low-speed 
streaks and one high-speed region form. In the case of two holes, three low- 
speed streaks and two high-speed regions form. Because of the particular 
direction of rotation of the counter-rotating vortices, only the middle low- 
speed streak will lift up, causing an unstable, inflectional velocity profile 
and bursting. Similarly, for three holes, two unstable streaks form, leading 
to two bursts. These arguments can be extended to an arbitrary number of 
spanwise suction holes. 

The artificial events resulting from withdrawing near-wall fluid in a 
laminar boundary layer through one, two, three or four suction holes were 
visualized using fluorescent dye seeping from the spanwise slot and a hori- 
zontal sheet of laser light at the surface of the plate. The results for a 
freestream speed of 20 cm/sec are shown in individual movie frames in 
Figures 6,7,8 and 9. The flow relative to the flat plate is from left to 
right, and the time from the onset: of suction is indicated on each movie 
frame. As expected, for a single suction hole, two stable streaks form and no 
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Figure 8, Four Streaks are Generated when Fluid is Withdrawn from Three 
Suction Holes. Only Two Low-Speed Streaks are Unstable, and 
Two Bursts are Generated 
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t = 1 .28 sec t = 1 .59 sec 

Figure 9. Withdrawing Fluid Through Four Suction Holes Results in the 


Generation of Five Low-Speed Streaks and Three Artificial Bursts 






bursting events are observed. Three low-speed streaks are marked by the bright 
dye accumulation in Figure 7 for the case of two suction holes, but only the 
middle one breaks into a burst. In Figures 8 and 9, two and three unstable 
streaks are observed, respectively, for the case of three and four suction 
holes . 

The hydrogen-bubble technique was also used to visualize the bursting 
events in a laminar boundary layer. Figure 10 shows a top view of the events 
resulting from withdrawing near-wall fluid through two suction holes separated 
in the spanwise direction. The hydrogen-bubble wire is perpendicular to the 
flow and parallel to the flat plate at a distance of 2 mm above the surface 
and 1 cm downstream of the suction holes. A horizontal sheet of laser light 
illuminates the bubbles. The low-speed regions are manifested as upstream 
kinks of the time lines, and longitudinal vorticity is indicated by the 
wrapping of these lines. The unstable region grows as it convects downstream 
and eventually breaks into a burst-like event. 

A side view of the events depicted in Figure 10 is visualized using a 
vertical hydrogen-bubble wire and a vertical sheet of laser light. The fluid 
is withdrawn from the laminar boundary layer through the two suction holes 
located 1 cm upstream of the hydrogen-bubble wire. The resulting events are 
shown in Figure 11 for six time intervals after the onset of the artificial 
event. The undisturbed Blasius profile can be reconstructed from the time 
lines shed from the hydrogen-bubble wire. The formation, lifting and 
stretching of the hairpin vortex that results from distorting the originally 
spanwise vortex lines is clearly indicated in the pictures shown in Figure 11. 
The vortex head reaches the outer region of the laminar boundary layer about 
1 sec after the suction has been initiated. The breakup of the vortex head is 
also indicated in the figure. 

Artificial bursts were also generated by pitching a flush-mounted delta 
wing as shown in the schematic in Figure 12. The negative angle of attack of 
30° leads to the formation of two longitudinal vortices. An unstable low-speed 
region is formed between the two vortices. The low-speed streak lifts up away 
from the wall, causing a strongly inflectional velocity profile, and rapidly 
breaks down into a completely random pattern much the same as the streaks 
generated by suction from spanwise holes. 
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An Artificial Burst Generated by Withdrawing Fluid from Two 

Suction Holes and Visualized Using Hydrogen Bubbles and 
a Horizontal Sheet of Laser Light 




Figure 1 1 . Side View of the Artificial Burst in a Laminar Boundary Layer 

Flow is Visualized Using Hydrogen Bubbles and a Vertical Sheet 
of Laser Light 
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Figure 12. Artificial Burst Generation. Resulting Streamlines When a 
Delta Wing, Originally Flush with the Surface, is Impulsively 
Pitched to Generate Negative Lift 
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4. Bursts in a Turbulent Boundary Layer 

A bursting event is, of course, a characteristic structure of a turbulent 
boundary layer. The artificial burst generators were used to generate burst- 
like events in a laminar flow as described in the previous section to allow 
easy visualization free of the random background that characterizes a turbulent 
flow. In this section, we describe the visualization of natural and artificial 
bursts in a turbulent boundary layer. The photographs in Figure 13 show the 
sequence of events leading to the breakup of a natural low-speed streak in a 
turbulent boundary layer. The wall region of the boundary layer is visualized 
by injecting fluorescent dye from the spanwise slot located just upstream of 
the view shown in the photographs. A horizontal sheet of laser at the wall 
clearly shows the natural low-speed streaks, manifested as bright regions of 
dye, with the typical spanwise spacing. In particular, the low-speed streak 
whose tail is marked with the vertical and horizontal arrows on the first 
photograph undergoes a large-amplitude oscillation at t = 0.28 sec (origin of 
time is arbitrary). At t = 0.41 sec, the breakup process has started, and it 
is completed at t = 0.81 sec, as shown in the last photograph in Figure 13. 

The breakup of a natural low-speed streak is compared to that of the 
artificially generated one in Figure 14. Here, the origin of time coincides 
with the onset of suction from two spanwise holes. In the first frame, the 
vertical and horizontal arrows indicate the position of the midpoint between 
the two suction holes. A low-speed streak is formed due to the pumping action 
of the longitudinal vortices that result from distorting the originally 
spanwise vortex lines. The artificially generated low-speed streak undergoes 
oscillation and breakup much the same as the natural one. The similarity 
between the natural and artificial events is striking, considering the fact 
that the natural phenomenon under consideration is random in space and time. 

A side view of the breakup of an artificially generated low-speed streak 
is shown in Figure 15. The turbulent boundary layer is visualized using 
fluorescent dye injected from the spanwise slot and a vertical sheet of laser 
light parallel to the flow. The dye accumulation near the wall indicates the 
formation of a low-speed region. In the first photograph in Figure 15, the 
streak lifts up as marked by the vertical arrow. The low-speed streak 
undergoes a large-amplitude oscillation and breaks up violently as shown in 
the sequence of movie frames. 
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Figure 13. Top View of a Turbulent Boundary Layer Visualized Using 
Fluorescent Dye and a Horizontal Sheet of Laser Light. 
Sequence of Events Leading to the Break-Up of a Natural 







t = 0,66 sec t = 0.81 sec 


Figure 14, Top View of a Turbulent Boundary Layer Visualized Using 
Fluorescent Dye and a Horizontal Sheet of Laser Light. 
Sequence of Events Leading to the Break-Up of an Artificially- 
Generated Low-Speed Streak 
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t = 0.08 sec 


FLOW 0 


t = 0.16 sec 


Figure 15. 


0.39 sec 


0.47 sec 


Side View of the Break-Up of an Artificially-Generated 
Low-Speed Streak. The Fluorescent Dye and the Vertical 
Sheet of Laser Light Show a Cut in the Hairpin Vortex 
Stretching and Lifting Away from the Surface 


28 



A different perspective of the natural and artificial bursts is obtained 
from side views using a vertical hydrogen-bubble wire and a vertical sheet of 
laser light. Figures 16 and 17 show the sequence of events leading to natural 
and artificially generated bursts, respectively. The formation of inflectional 
velocity profiles is clearly indicated in both figures. These profiles are 
inviscidly unstable, and the flow rapidly breaks down into a completely random 
pattern. It is interesting to note the different time scales associated with 
the bursting event in laminar and turbulent boundary layers (e.g., Figures 11 
and 17). This is due to the different friction velocities and thus different 
viscous scales in the two boundary layers. 
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t = 0.03 sec 


FLOW 0 


0.1 1 sec 



Side View of a Natural Burst in a Turbulent Boundary Layer 
Using Hydrogen Bubbles and a Vertical Sheet of Laser Light 
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Figure 17. Side View of an Artificial Burst Generated by Withdrawing 
Fluid from the Turbulent Boundary Layer Using Two Suction 
Holes. The Flow is Visualized Using Hydrogen Bubbles and 
a Vertical Sheet of Laser Light 
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5. Pattern Recognition of Bursts and Streaks 

The visualization results presented in Sections 3.3 and 3.4 indicate the 
sequence of events associated with the generation of an artificial burst in 
laminar and turbulent boundary layers, respectively. The observed flow 
patterns agree qualitatively with observations of naturally occurring bursts 
and streaks in a turbulent boundary layer. The length and time scales of the 
artificial events also agree with those for natural structures. For quantita- 
tive comparison between natural and artificial bursts, hot-film probes were 
used together with pattern recognition algorithms aimed at detecting low-speed 
regions and bursting events. 

The miniature hot-film probes were first used to obtain mean and rms 
strearawise velocity profiles. In the laminar case, the mean velocity profiles 
were found to agree with the Blasius solution for a zero-pressure-gradient 
boundary layer. Similarly, the mean and rms velocity profiles in the turbulent 
case agreed well with those for a standard flat plate boundary layer. The 
mean profiles always had well-defined linear, logarithmic and wake regions. 
Clauser's (1956) method was used to determine the friction velocity used for 
normalizing the data and in the pattern recognition algorithms described below. 

The hot-film probes used to detect the bursts and streaks were placed 100 
to 1000 wall units downstream of the artificial burst generator, and at a 
height above the wall in the range of 10 to 20 wall units (a wall unit was 
about 0.1 mm at a freestream speed of 20 cm/sec and a streamwise distance from 
the leading edge of 80 cm). The instantaneous longitudinal velocity signal 
was mapped at these locations. In the laminar case, the initiation of an event 
using either the suction holes or the pitching delta wing was always associated 
with the passage of a velocity "spike" before the flow broke down into turbu- 
lence. These velocity spikes are consistent with the passage of the head of a 
hairpin vortex and are similar to the spikes observed by Klebanoff et al . 

(1962) in their classical vibrating ribbon experiment. In the case of a turbu- 
lent boundary layer, the artificial event caused a large streamwise velocity 
acceleration in the instantaneous velocity records similar to the accelerations 
observed by Blackwelder & Kaplan (1976) for natural bursts in a turbulent 
boundary layer. For a more objective assessment, the VITA technique was used 
to detect the passage of a burst, and the streak detection algorithm was used 
to recognize the formation and the evolution of a low-speed region. Both 
pattern recognition algorithms are reproduced in Appendix I. 
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The two pattern recognition algorithms, which have been described briefly 
in Section 3.2, search for recognizable patterns or structures. The burst 
detector searches for periods of high acceleration without necessarily larger 
fluctuations than the background turbulence. As shown by Blackwelder & Kaplan 
(1976), such periods are associated with events characterized by a high degree 
of coherence in time and in the direction normal to the wall and having a con- 
ditionally averaged Reynolds shear stress that is an order of magnitude greater 
than its conventionally averaged value. The detection scheme is then based on 
the dynamic properties of the organized structure. Similarly, the streak 
detector searches for low-speed regions near the wall relative to the 
background mean velocity. 

Whenever the artificial burst generator was triggered, a low-speed streak 
and a burst were detected, after an appropriate time delay, by the two pattern 
recognition algorithms. This confirms our hypothesis that the artificial 
event is dynamically similar to a natural burst. 

An important quantity to further confirm the dynamic similarity between 
natural and artificial events is the conditionally averaged Reynolds stress 
<uv>. This statistical quantity was not directly measured in the present pre- 
liminary investigation, but will be measured during the second phase of this 
work. 
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6 . Summary 

The objective of the Phase I work was to demonstrate the feasibility of 
generating artificial bursts in laminar and turbulent boundary layers. Experi- 
ments were conducted in FLOW's 18-meter towing tank. Zero-pressure-gradient 
laminar or turbulent boundary layers were generated by towing a flat plate in 
the water channel. Artificial bursts were generated by either withdrawing 
near-wall fluid from one, two, three or four minute holes separated in the 
spanwise direction, or by pitching a miniature delta wing that is flush-mounted 
to the wall of the plate. The resulting events were visualized using fluores- 
cent dyes or hydrogen bubbles illuminated by a sheet of laser light projected 
in the desired plane. Both the suction and the pitching of the delta wing 
generated streamwise vorticity and one or more low-speed streaks. The unstable 
streaks lifted away from the wall, causing an inflectional velocity profile, 
oscillated and broke up into a burst. 

Hot-film probes were used to measure the instantaneous longitudinal 
velocity, and the resulting signals were analyzed using two pattern recognition 
techniques that detect bursts and low-speed streaks. The two detectors consis- 
tently recognized a burst and a low-speed streak whenever the artificial burst 
generator was triggered. 

The artificial events generated in this study are similar qualitatively 
and quantitatively to natural bursts in a turbulent boundary layer. We are 
now in a position to examine in detail the artificial bursts that occur at a 
given streamwise location in the boundary layer and at controlled times. The 
topological properties of the artificial bursts will be educed from phase- 
locked hot-film measurements. Means of manipulating natural and artificial 
bursts for the purpose of boundary layer control and drag reduction will be 
contemplated during the second phase of this investigation. 

By all available measures, Phase I of this project should be considered 
highly successful. We have been able to establish the proof-of-concept : the 

two-hole suction and the pitching delta wing do indeed produce longitudinal 
vortices and the associated longitudinal streaks of high-speed and low-speed 
fluids separated in the spanwise direction; the low-speed streak becomes 
unstable and breaks down in a way which appears to be identical with the 
bursting process. But Phase I has yielded a lot more: the dye injection 

method is successful in revealing the artificial bursting phenomenon; the 
hydrogen-bubble visualization is particularly helpful in obtaining a clear 
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perception of the time-evolution of the three-dimensional event; a series of 
holes separated in the spanwise direction can produce a series of simultaneous 
parallel streaks and bursting, etc., as discussed in Section 3. These results 
are more than what could be realistically expected during Phase I and are 
highly encouraging. They indicate that far greater successes are likely in 
Phase II. 

The results of the Phase I work will be presented at the 38th Annual 
Meeting of the American Physical Society, 24-26 November, 1985, at the AIAA 
24th Aerospace Sciences Meeting, 6-8 January 1986, and at the IUTAM Symposium 
on Fluid Mechanics in the Spirit of G. I. Taylor, 24-28 March 1986. A 
manuscript describing our experimental results will be sent to the Journal of 
Fluid Mechanics. 
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Appendix I: Pattern Recognition Algorithms 
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n n n tn -fc> cj n non 


c #**iHHHHt*** PROGRAM TO COMPUTE THE CONVENTIONAL STATISTICS 


C *********** PROGRAM TO COMPUTE THE CONVENTIONAL STATISTICS 
C AND STORE THE RESULTS IN AN INFORMATION BLOCK 

C 

C MODIFIED BY PAT MCCAFFERTY B50530 TO MAKE ,, IBA M INTEGERS TO 
C PREVENT RMS BECOMING NEGATIVE WHEN T EXCEEDS 30 SECONDS 
C WITH DT = . 001 (32000 PTS INTEGER SHORT LIMIT) 

C 

PARAMETER (LEN=11> 

REAL A ( 32 ) , UBAR(LEN), RMS(LEN) , C(4, LEN), B(128) 

REAL*B SUMU(LEN), SUMSQ(LEN) 

REAL*B AA, BB 
CHARACTER -*32 FNAME 
CHARACTER *5 RUNO 
INTEGER*4 KCHN ( 5 ) , LCHN ( 32 > , I BA 
INTEGER*2 IU(4096> 

DATA SUMU/LEN*0. /, SUMSQ/LEN*0. /, UBAR/LEN*0. /, RMS/LEN*0. / 

********** READ INDX FILE 

PRINT *, 'ENTER LAST FIVE CHARACTERS OF RUN NAME' 

READ(1, ' ( A5 ) ' ) RUNO 
FNAME= ' INDX/COAT ' //RUNO 
PRINT *, ' FILE NAME IS '.FNAME 

OPEN ( 40, FILE=FNAME, STATUS* 'OLD ', RECL=256, ERR=5. 

*FORM= 'UNFORMATTED ', ACCESS* 'DIRECT ' ) 

READ(40, ERR = 5) ( IU < K ) , K=1 , 256 ) 

CLOSE (40) 

KCHN ( 1 ) = IU(181) / *TOTAL NUMBER OF CHANNELS RECORDED 

IF ( KCHN ( 1 ) . GT. LEN ) STOP 'PROGRAM ONLY ACCEPTS 11 CHANNELS' 
KCHN ( 2 ) = IU( 182) /*1ST CHN OF ADC'S USED 

CALL RN2P( IU(9), DT, INTL( 1 ) ) 

WRITE (1,3) DT 

FORMAT ( /, ' DIGITIZING INTERVAL =',FB. 5, ' SECONDS. ',/) 

A ( 21 ) = DT /^DIGITIZING INTERVAL 

CALL RN2P< IU( 159), UOO, INTL(l)) 

WR ITE (1.4) UOO 

FORMAT ( ' RECORDED VALUE OF UOO =',F6. 1, ' CM/SEC. ',/) 

A (22) = UOO 
GO TO 10 

PRINT ■*, ' ERROR ENCOUNTERED IN READ ING ', FNAME 

******** read THE COEFFICIENT FILE 

10 FNAME* 'COEF/COAT '//RUNO 

C PRINT *, ' COEFFICIENT FILE IS '.FNAME 

0PEN(31 , FORM* 'UNFORMATTED ', STATUS* 'OLD ' , FILE=FNAME ) 

READ ( 31 ) B 

CALL RN2P(B, B, INTL( 128) ) 

CLOSE (31 ) 

IC=4*KCHN(2> /*BIAS COEF TO FIRST CHAN RECORDED 
DO 20 K = 1, KCHN ( 1 ) 

LCHN(K) =K+KCHN(2)-1 
DO 20 J = 1,4 
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C *********** PROGRAM TO COMPUTE THE CONVENTIONAL STATISTICS 


IC = IC + 1 
20 C< J, K) = B( IC) 

PRINT *, ' CHANNEL A B C D' 

WR I TE < 1 > 23 ) (LCHN(L) . ( C ( J, L ) , J=1 , 4 ) , L=l, KCHN< 1 ) ) 

23 FORMAT ( I5< 6X. 4G10. 3) 

******** INPUT CHANNEL NUMBER OF HOT-FILMS 

PRINT *, 'TOTAL NUMBER OF CHANNELS RECORDED = ',KCHNd) 

CALL TNOUA( ' ENTER CHANNEL NUMBER OF FIRST HOT-FILM: ',42) 

READ(1,*) KCHN ( 3 ) 

CALL TNOUA( ' ENTER TOTAL NUMBER OF HOT-FILMS RECORDED: ',44) 

READ(1,*) KCHN ( 4 ) 

I FILM = KCHN ( 4 ) 

********** SELECT THE TIME INTERVAL TO ANALYSE 

RECL = INT (4096/KCHNd ) )*DT /*RECORD LENGTH IN SECONDS 
CALL TNOUA ( ' ENTER THE STARTING T I ME C SEC 3 : ',32) 

READd,*) TSTART 
ISTART = TSTART/RECL+1 
TSTART =( ISTART-1 )*RECL +DT 
PRINT *, 'TSTART =', TSTART 

CALL TNOUA ( ' ENTER THE STOPPING TIMECSEC3: ',32) 

READd,*) TSTOP 
ISTOP = TSTOP/RECL+1 
TSTOP =ISTOP * RECL 
PRINT *, 'TSTOP =', TSTOP 

*********** COMPUTE THE MEAN AND RMS VALUES OF ALL CHANNELS 

FNAME= 'DATA/COAT '//RUNO 
ICH = KCHN ( 1 ) 

C2V=1. /409. 6 
C2V2=C2V*C2V 
C2V3=C2V2*C2V 
JBLK = 0 
IBA = 0 

LRCL=( 4096/ ICH)* ICH 

OPEN ( 29, FORM= 'UNFORMATTED', STATUS 3 'OLD ' , FILE=FNAME, ERR=901) 

40 READ ( 29, END=1 10. ERR=902 > IU 
JBLK = JBLK + 1 
IF( JBLK. LT. ISTART) GO TO 40 
IF ( JBLK. GT. ISTOP) GO TO 40 

C WRITE! 1,11) (IU(M), M«l, LEN) 

11 FORMAT ( 519) 

C ADD THE PRESENT DATA TO THE SUMMED STATISTICS 

IBA = IBA + 1 

DO 55 K=KCHN(3), KCHN ( 3 ) +KCHN ( 4 ) — 1 
AA=0. DO 
BB=0. DO 
C0F1=C<1, K) 

C0F2=C < 2, K ) *C2V 
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c *-*•#**■»***** PROGRAM TO COMPUTE THE CONVENTIONAL STATISTICS 


C0F3=C <3, K)*C2V2 
C0F4=C<4, K)*C2V3 
DO 50 J=K, LRCL, ICH 
V=FLOAT < IU(U) ) 

U=C0F1+V* ( C0F2+V-* < C0F3+V*C0F4 ) ) 

AA=AA+U 

50 BB=BB+U*U 

SUMU ( K ) =SUMU ( K ) +AA 
55 SUMSQ ( K ) =SUMSQ ( K ) +BB 

C PRINT *, ' FINISHING ', JBLK, ' PASS' 

GO TO 40 
C 

C -»******•# NORMALIZE THE DATA 
C 

110 CLOSE ( 29 ) 

NORM = 4096/ ICH 

XNORM = IBA*NORM /-^NORMALIZATION FACTOR 

DO 200 K = KCHN<3). KCHN ( 3 ) +KCHN ( 4 ) -1 

UBAR < K ) = SUMU ( K ) /XNORM 

XX = SUMSQ<K) /XNORM 

IF< XX - UBAR(K)*-»2) 190,195,195 

190 PRINT *, ' NEGATIVE SORT ON THE ', LCHN( K ) , ' CHANNEL ==> RMS=0. ' 
RMS(K) = 0. 

GO TO 200 

195 RMS(K) = DSQRT ( SUMSQ (K ) /XNORM - UBAR(K)**2) 

200 CONTINUE 
C 

C ■*■*****■»* OUTPUT THE DATA 
C 

A ( 23 ) = JBLK 
A ( 24 ) = ISTART 
A ( 25 ) = ISTOP 

0PEN(32, FILE= 'BURST. DATA '//RUNO, FORM= 'UNFORMATTED ' ) 

WRITE(32)A, C, UBAR, RMS, KCHN, FNAME /*WRITE INFORMATION ONTO DISK 
CLOSE (32) 

TREC = RECL*IBA 
XREC = RECL* JBLK 
WRITE( 1,209) TREC, XREC 

209 FORMAT ( /, ' THE MEAN STATISTICS F0R',F6. 1, ' SEC OF DATA ARE: ' 

*, /, 5X'(THE TOTAL AVAILABLE RECORD IS',F6. 1, ' SEC. )') 

WRITE! 1,211) <LCHN< J), UBAR (J), RMS (J), J=l, ICH) 

211 F0RMAT(/,3X, 'CHANNEL # UBAR ', 9X, 'RMS ', /» 12(4X, 17, 2F12. 2. / ) ) 




CALL 

EXIT 



C 

901 

DATA 

FILE 

OPEN 

ERROR 

C 

902 

DATA 

FILE 

READ 

ERROR 


901 

STOP 

901 




902 

STOP 

END 

902 
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C 


PROGRAM (STREAK) 


C PROGRAM (STREAK) 

C ******** COMPUTES THE LOCATIONS OF LOW SPEED REGIONS FROM THE 3 
C SPANWIZE HOT-FILMS USING THE SECOND DERIVATIVE OF 

C DU/DZ OR A THRESHOLD ON THE VELOCITY DATA AND 

C PLOTS THE HISTROGRAM OF THE STREAK LENGTHS. 

C THIS PROGRAM ASSUMES THE MEAN AND RMS STATISTICS 

C HAVE BEEN CALCULATED AND STORED ON THE DISK. 

C > CHANGING 'LENH' VARIES THE ABSCISSA OF THE HISTOGRAM 

C > ' XAM ' CONTROLS THE NOMALIZATION OF THE HISTOGRAM ORDINATE. 

C 

PARAMETER ( LEN=1 1 ) 

PARAMETER (LENH=60> 

REAL A (32), UBAR(LEN), RMS(LEN), C(4, LEN), B ( 128) , OLDT, NOWT 
REAL UE (4096 ) , U(3), UT (3) 

INTEGER*2 IU(4096), HIST (LENH), ICH(3) 

INTEGER*4 I TOT, ISAV, KCHN ( 5 ) 

CHARACTER AST, BLNK, LINE(60) 

CHARACTER*32 FNAME 
CHARACTER*5 RUNO 

DATA AST/ '*'/, BLNK/ ' '/, LINE/60* HIST/LENH*0/ 

******** READ IN STORED INFORMATION FROM THE DISK 

PRINT *, ' ENTER LAST FIVE CHAR OF RUN NAME : ' 

READ ( 1 , ' (A5> ' ) RUNO 

OPEN < 32, FILE= 'BURST. DATA '//RUNO, FORM = 'UNFORMATTED ' , 

* STATUS= 'OLD ' ) 

READ ( 32 ) A, C, UBAR, RMS, KCHN, FNAME /* INPUT STORED INFORMATION 

CLOSE (32) 

******** SELECT THE 3 CHANNELS FOR THE STREAK ANALYSIS 
IF ( KCHN ( 4 ) . GT. 2 ) GO TO 5 

PRINT *, ' LESS THAN 3 CHANNELS OF H. F. RECORDED !!!' 

STOP 

WR ITE (1,6) KCHN ( 3 ) 

FORMAT ( ' THE FIRST H. F. SIGNAL IS ON CHANNEL #',I3, '. ') 

CALL TNOUA ( ' ENTER CHANNEL # OF FIRST H. F. TO ANALYSE : ',43) 

READ(1,*) ICH(l) 

ICH ( 2 > = ICH(l) + 1 /* ASSUMES H. F. CHANNELS ARE SEQUENTIAL 

ICH ( 3 ) = ICH ( 1 ) +2 /» ASSUMES H. F. CHANNELS ARE SEQUENTIAL 

PRINT *, ' THE STREAK DETECTION WILL USE CHANNELS : ICH 

• 

******** SET UP PARAMETERS FOR DETECTING THE STREAKS 

NCH = KCHN ( 1 ) /* TOTAL # OF CHANNELS RECORDED 

RECL= INT ( 4096 /NCH ) *A ( 21 ) /* A(21)=DT 

THRES =1. /* THRESHOLD USED FOR VELOCITY 

DO 10 K = 1,3 

10 ICH(K) = I CH ( K ) -KCHN ( 2 ) +1 /* BIAS ICH TO 1ST CHANNEL INDEX 

UT ( 1 ) = UBAR ( ICH ( 1 ) ) - 0. 5*THRES*RMS( ICH( 1 ) ) /* VELOCITY THRESHOLD 

UT ( 2 ) = UBAR ( ICH ( 2 ) ) - 1 . 0*THRES*RMS ( ICH ( 2 ) ) /* VELOCITY THRESHOLD 

UT ( 3 ) = UBAR ( ICH ( 3 ) ) - 0. 5*THRES*RMS ( ICH ( 3 ) ) /* VELOCITY THRESwni n 
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c 


PROGRAM (STREAK) 


ITOT = 0 

/* 

OLDT = 0. 

/* 

NOWT = 0. 

/* 

GAMMA = 0. 



TIME COUNTER 

STATUS WORD OF PREVIOUS DATA POINT 
STATUS WORD OF PRESENT DATA POINT 


IT = 0 
DT = A ( 21 ) 

UOO = A (22) 

NBLK = A (23) 

ISTART = A ( 24 ) 

ISTOP = A ( 25 ) 

CALL TNOUA ( ' ENTER 
READ(1,*) UTAU 

PRINT *, 'STREAK STARTING TIME 


/* 

/* 

/* 

/* 

/* 

/* 


COUNTS NUMBER OF STREAKS 
DIGITIZING INTERVAL 
TOWING SPEED 

TOTAL # OF C4096I BLOCKS OF DATA 
ISTART=FIRST BLOCK OF DATA TO ANALYSE 
ISTOP=LAST BLOCK OF DATA TO ANALYSE 


THE VALUE OF UTAU [CM/SEC 3 


34) 


STREAK STOPPING TIME' 


C 

C ******** BEGIN MAIN LOOP 
C 

IBLK = NBLK 

LBLN = ( 4096/NCH ) *NCH 

OPEN ( 29/ FORM= 'UNFORMATTED ', STATUS= 'OLD ' , FILE=FNAME ) 

DO 210 IB = 1, IBLK 

READ ( 29, ERR =901 ) IU 

IF( IB. LT. ISTART) GO TO 210 

IF( IB. GT. ISTOP) GO TO 210 

CALL LINEAR3 ( IU, UB, C, ICH( 1 ) , NCH > 

DO 200 J = ICH( 1 >, LBLN, NCH 
ITOT = ITOT + 1 
DO 30 K = 1, 3 
30 U(K) = UB ( J+K-l ) 

CR = U ( 3 > - 2. *U ( 2 ) + U(l) /* C = CURVATURE = 2ND DERIVATIVE 
IF(CR) 35,35,40 /* IF C>0 ==>MINIMUM==>POSSIBLY STREAK 

35 I F ( U ( 1 ) - UT ( 1 ) ) 3B, 39,150 /* STRONG DEFECT==> STREAK 

38 IF ( U ( 3 ) - UT (3 ) ) 40,40,150 /* STRONG DEFECT==>STREAK 

40 IF(U(2) - UT ( 2 ) ) 60,60,150 /* STRONG DEFECT==>STREAK 

C 

C ******** STATEMENT 60==> THERE IS A SLOW SPEED REGION 
C 

60 NOWT = 1. 

GAMMA = GAMMA + 1. 

IF(OLDT. EQ. NOWT) GO TO 200 
C ******** THIS IS A NEW LOW SPEED REGION 


DUMIB=IB 
DUM J= J 

Tl=( (DUMIB-1. ) + ( ( DUM J— 1 . ) /4096. ))*RECL 
IT = IT + 1 
ISAV = ITOT 
GO TO 190 


C 

C ******** STATEMENT 150==> THERE IS NO LOW SPEED REGION 
C 

150 NOWT = O. 

IF(OLDT. EQ. NOWT) GO TO 200 
C ******** THIS IS THE END OF THE LOW SPEED REGION 
DUM I B= IE 
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c 


PROGRAM (STREAK) 


DUM J= J 

T2=< (DUMIB-1. ) + ( ( DUM J- 1 . )/4096. >)*RECL 
WRITEd, '(2<5X, F8. 4) > ') Tl, T2 
LPOS = ITOT - ISAV 
IF(LPOS. GT. LENH) LPOS = LENH 
HIST (LPOS) = HIST ( LPOS ) + 1 
190 OLDT = NOWT 
200 CONTINUE 
210 CONTINUE 
CL0SE(29) 

C 

C ******** COMPUTE THE MEAN STATISTICS 

TREC = ITOT*DT /* TOTAL RECORD LENGTH 

PLS = UTAU*UTAU/0. 01 /* NORMALIZING CONSTANT 

IF( IT . EG. 0) WRITEd, 220) 

IF( IT . EQ. O) GO TO 901 

220 FORMAT ( ' **** NO STREAKS WERE DETECTED **** ') 

ALNGTH = GAMMA*DT/IT*PLS /* AVERAGE LENGTH IN T+ 

GAMMA = GAMMA /I TOT* 100. /* INTERMITTENCY FACTOR 

FREQ = IT/TREC /* FREQUENCY 

FPLS = FREQ/PLS /* F+ 

WRITEd. 251 ) RUNO, IT, TREC, GAMMA, FREQ, FPLS, ALNGTH 
251 FORMAT(///, 10X, 'FOR RUN # '.AS?', ',15, ' STREAKS', 

*' WERE FOUND IN A',/, 15X, 'TOTAL RECORD 0F',F7. 2, 

* 'SECONDS. ', /, 20X, 'GAMMA =',F8. 2, 20X 

*, 'FREQUENCY =', F7. 1, ' HZ', /, 20X, 'F+', SX, '=', F8. 4, / 

*, 20X, 'LENGTH T+ =', F8. 1, /) 

C 

C ******** FIND THE MAXIMUM IN THE HISTOGRAM 
C 

MAX = 0 

DO 300 I = 1 , LENH 
IF(HIST(I). LT. MAX) GO TO 300 
MAX = HIST ( I ) 

300 CONTINUE 

340 XAM = FLOAT ( MAX ) /60. 

C XAM = 1. /* XAM = 1. ==>UNNORMAL I ZED HISTOGRAM 

C 

C ******** PLOT THE HISTOGRAM 
C 

DO 520 I = 1, 60 
520 LINE ( I ) = '-' 

WRITEd. 551) LINE /* PRINT THE ORDINATE AXIS 

551 FORMAT (//, 11X,6(9X, ': '),/, ' LEN # ',60A1) 

DO 560 J = 1. 60 

560 LINE(U) =* BLNK /* ZERO THE LINE ARRAY 

DO 600 I = 1 , LENH 

STL = I*DT*PLS /* STREAK LENGTH IN T+ 

UH = HIST ( I ) /XAM +0.5 /* JH = SCALED ORDINATE OF HISTOGRAM 

IF( JH. NE. O) GO TO 575 

WRITEd, 566) STL, HIST(I) /* OUTPUT HISTOGRAM IF UH=0 

566 FORMAT ( F5. 1. 14. ' : ' ) 

GO TO 600 
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non 


C 


PROGRAM (STREAK) 


575 

IF( JH. GT. 60) JH = 60 
DO 580 J = 1, JH 



580 

LINE ( J ) = AST 

/* 

PLACE INTO LINE ARRAY 


WRITE( 1, 581 ) STL, HIST( I ), LINE 

/* 

OUTPUT HISTOGRAM 

581 

F0RMAT(F5. 1, 14, ' :',60A1> 

DO 590 J = 1, JH 



590 

LINE ( J ) = BLNK 

/* 

ZERO THE LINE ARRAY 

600 

CONTINUE 



901 

STOP 901 




END 

C 

SUBROUTINE LINEAR3(IU. U, C,LCH, NCH ) 

**«##*** LINEARIZES 3 DATA CHANNELS OUT OF A TOTAL NCH CHANNELS 
IU = INPUT INTEGER ARRAY U = OUTPUT REAL ARRAY 
C = COEFFICIENT ARRAY LCH = STARTING CHANNEL 
INTEGERS IU ( 4096 ) 

REAL U ( 4096 ) » C(4, 11 ) 

LCH3 = LCH + 2 
DO 100 L = LCH, LCH3 
DO 100 I = L, 4096, NCH 
X = IU(I-1 )/409. 6 
100 U(I) = C ( 1 , L ) + X* ( C ( 2, L ) + X* ( C ( 3, L ) + X*C(4,L)>) 

RETURN 

END 
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C ********** THIS PROGRAM COMPUTES THE BURSTING TIMES AND 


********** THIS PROGRAM COMPUTES THE BURSTING TIMES AND 
REJECTS THE DECELERATIONS. 

PARAMETER ( LEN=1 1 ) 

REAL A ( 32 ) > UBAR(LEN), RMS(LEN), OLDT, NOWT, C<4, LEN > , UB(5096) 
INTEGER*4 ITOT, ISAVE, KCHN < 5 ) 

INTEGER*2 IU< 5096 > , HIST< 40 ) 

REAL*S VITAE1 , VITASQ, U, Ul. AA, BB, V. W 

CHARACTER*32 FNAME 

CHARACTER *5 RUNG 

COMMON /DATA/ IU 

DATA HIST/40*0/ 

********* READ IN STORED INFORMATION FORM DISK 

PRINT *, 'ENTER LAST FIVE CHAR OF RUN NAME : ' 

READ ( 1. ' <A5> ' ) RUNO 

OPEN( 32< FILE= 'BURST. DATA ' //RUNO. FORM= 'UNFORMATTED ' , 

* STATUS= 'OLD ' ) 

READ ( 32 ) A, C, UBAR, RMS, KCHN, FNAME /*INPUT STORED INFORMATION 
CLOSE < 32) 

******* SELECT THE TRIGGER CHANNEL FOR THE VITA TECHNIQUE 
IF ( KCHN ( 4 ) — 1 ) 5,9,7 

CALL TNOUA( ' ENTER CHANNEL NUMBER OF FIRST H. F. : ',38) 

READd,*) KCHN ( 3 ) 

CALL TNOUA( ' ENTER CHANNEL # OF H. F. TO ANALYSE: ',38) 

READd,*) ICH 
GO TO 10 

ICH = KCHN ( 3 ) /*ICH = CHANNEL # OF H. F. 

| PRINT *, ' CHANNEL #', ICH, ' WILL BE ANALYSED. ' 

********* SET UP PARAMETERS FOR VITA TECHNIQUE 

DT = A ( 21 ) /*DIGITIZING INTERVAL 

UOO = A ( 22 ) / *TOWING SPEED 

NBLK = A ( 23 ) /*TOTAL NUMBER OF C40961BL0CKS OF DATA 

ISTART = A ( 24 ) /*ISTART = FIRST BLOCK TO ANALYSE 

I STOP = A < 25 > /*ISTOP = LAST BLOCK TO ANALYSE 

NCH = KCHN ( 1 ) 

RECL=INT ( 4096/NCH)*DT 

ICH= ICH-KCHN ( 2 > + 1 /*BIAS ICH TO 1ST CHAN INDEX 

AVGU = UBAR (ICH) 

URMS = RMS < ICH) 

PRINT *, 'UOO = ',U00, 'UBAR =',AVGU. 'RMS = ',URMS 
CALL TNOUA ( ' ENTER THE VALUE OF UTAUCCM/SEC I : ',35) 

READd,*) UTAU 

DELT = 0. 01/UTAU/UTAU 

NSAMP = 10. *DELT/DT 

NSAMP = ( NSAMP /2) *2 /*MAKE NSAMP EVEN 
PRINT *, ' NSAMP = ', NSAMP 
NS = NSAMP /3*NCH 
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C ********** THIS PROGRAM COMPUTES THE BURSTING TIMES AND 


C 


FNSMP = 1. /NSAMP 
THRES =1.0 
ULEV = THRES*URMS**2 
GAM = 0. 

OLDT = O. 

TRANS = 0. 

ITOT = 0 
VITAE1 = 0. DO 
VITASQ = 0. DO 


/♦TRIGGER LEVEL -1.2 WAS USED BY RFB AND REK 
/♦ZERO INTERMITTENCY 

/♦OLDT IS THE FLAG FOR THE LAST DATA POINT 
/♦# OF TRANSITIQNS( I. E. BURSTS) 

/♦ITOT COUNTS THE RECORD LENGTH 
/♦VITA SUM FOR MEAN VALUE 
/♦VITA SUM FOR MEAN SQUARE VALUE 


PRINT *, 'BURST STARTING TIME BURST STOPPING TIME' 

******** BEGIN THE MAIN LOOP 


IBLK = NBLK 

LBLN= < 4096 /NCH ) ♦NCH+ 1 000 

OPEN ( 29, FORM= 'UNFORMATTED '. STATUS* 'OLD ' , FILE=FNAME ) 

DO 210 IB = 1, IBLK 

READ (29, ERR=901 ) ( IU< I >, 1 = 1001, 5096) 

IF( IB. LT. ISTART) GO TO 210 

IF( IB. GT. ISTOP ) GO TO 212 

CALL LINEAR ( IU( 1001 ), UB( 1001 ), C, ICH, NCH ) 

DO 200 IC = ICH+1000, LBLN, NCH 

ITOT = ITOT + 1 /♦INCREMENT POSTITON COUNTER 

U = UB ( IC ) -AVGU 

U1=U 

VITAE1 = VITAE1 + U /♦SHORT TIME AVERAGE OF U 

VITASQ = VITASQ + U*U /*SHORT TIME AVERAGE OF U^*2 

IF ( ITOT. LE. NSAMP) GO TO 200 
JC = IC - NSAMP*NCH 
U = UB ( JC ) —AVGU 

VITAE1 = VITAE1 - U /*SUBTR ACT OLD DATA POINT 

VITASQ = VITASQ - U*U /♦SUBTRACT OLD DATA POINT 
23 VARNCE = VITASQ*FNSMP - ( VI TAE1 ♦FNSMP )«*2 /♦COMPUTE VARANCE 

IF ( VARNCE + .0001) 100,110,110 

100 WRITEC 1,101) VARNCE, ITOT, IC, U, VITAE1, Ul, VITASQ 

101 F0RMAT(F10. 5, 218, 4F12. 5, /) 

XSTOP = XSTOP + 1 

IFCXSTOP.GT. (NSAMP+10. )) STOP 'VARIANCE WAS NEGATIVE' 

110 CONTINUE 

IF( VARNCE - ULEV) 150, 160, 160 
C ♦*♦♦*♦♦♦* STATEMENT #150 IS FOR I(T) =0 
150 NOWT = 0. 

IF ( NOWT - OLDT) 155, 190, 190 

C STATEMENT 155 IS EXECUTED WHEN EXITING FROM A BURST 

155 TRANS = TRANS + 1 /♦COUNT THIS LAST BURST 

DUMIB=IB 
DUMIC=IC— 1000 

T2=< (DUMIB-1 . )+( (DUMIC-1. )/4096. ))^RECL 
WRITEC1. '<2<5X, F0. 4) ) ') T1,T2 
I SAVE = O 
UMIN = 1000. 

UMAX = -1000. 

DO 157 I = JC, IC, NCH 


48 



c ********** THIS PROGRAM COMPUTES THE BURSTING TIMES AND 


IF(UB( I ). LT. UMIN) UMIN = UB ( I ) 

IF(UB( I >. GT. UMAX) UMAX = UB ( I ) 

157 CONTINUE 

IH » UMAX - UMIN +1.5 
IF( IH. LT.l) IH = 1 
IF< IH. GT. 40) IH = 40 
HIST(IH) = HIST ( IH) + 1 
GO TO 190 

C ******** IS THIS AN ACCELERATION OR DECELERATION? 

160 IFCUB(IC-NS) — UB ( JC+NS ) ) 150,150.170 /^ELIMINATE DECELERATIONS 
C *********** STATEMENT #170 IS FOR I(T) = 1, I. E. INSIDE A BURST 
170 NOWT = 1. 

GAM = GAM + 1. 

IF(NOWT - OLDT) 190, 190, 180 

180 ISAVE = ITOT /+SAVE BEGINNING COUNTER OF THIS BURST 

DUMIB=IB 
DUMIC=IC-1000 

Tl = ( (DUMIB-1. > + ( ( DUMIC-1 . >/4096. ))*RECL 
190 OLDT = NOWT 
200 CONTINUE 

C ********** MOVE REMAINING DATA TO THE BEGINNING OF THE BUFFER 
LSPT = ( 4096 /NCH ) +NCH + 1 
JC = 0 

DO 205 IC = LSPT, LSPT+1000 
JC = UC + 1 
205 UB( JC) = UB( IC) 

210 CONTINUE 

212 TREC = (ITOT - NSAMP ) *DT /*TOTAL RECORD LENGTH 

BURSTS = TRANS /^NUMBER OF BURSTS 

FREQ = BURSTS/TREC 

FIN = FREG*DELT 
GAMMA = GAM/ 1 TOT* 100. 

WRITE! 1, 213) THRES, BURSTS, GAMMA, FREQ, FIN 

213 FORMAT!/, ' THE 3URSTING STATISTICS FOR THRESHOLD = ', 

1F4. 1, ' ARE: ', /, 10X, 'NUMBER OF BURSTS =', F10. 0, 

2/, 10X, 'INTERMITTENCY FACTOR =',F8. 2, ' 7.' 

3/, 10X, 'FREQUENCY =', F19. 2, ' HZ ', /, 10X, 'FREQ+ =', F26. 5, /) 

CALL HSTER (HIST, 40) 

CLOSE ( 29 ) 

CALL EXIT 

C 901 DATA READ ERR 
901 STOP 901 
END 
C 

SUBROUTINE LINEAR! IU, U, C, ICH, NCH) 

INTEGER*2 XU (4096) 

REAL U( 4096 ) , C (4, 11 ) 

DO 100 I = ICH, 4096, NCH 
X = IU! I )/409. 6 

100 U(I) = C( 1, ICH)+X*(C(2, ICH)+X*(C(3, ICH)+X*C(4, ICH) ) ) 

RETURN 

END 

C 
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c ********** THIS PROGRAM COMPUTES THE BURSTING TIMES AND 


SUBROUTINE HSTER (HIST, N) 

INTEGER*2 HIST(N) 

CHARACTER LINE(40), BLNK, AST 
DATA BLNK/' '/.AST/'*'/ 

C 

C LOCATE FIRST AND LAST POINT IN HISTOGRAM 

C 

IMIN = N 
DO 20 I = 1, N 
IFCHIST ( I ) . LE. 0 ) GO TO 20 
IMIN = I 
GO TO 25 
20 CONTINUE 

25 IMAX = 0 

DO 30 I = 1, N 

IF(HIST(N+1-I). LE 0) GO TO 30 
IMAX = N + 2 - I 
GO TO 35 
30 CONTINUE 

35 WRITE( 1, 39) 

39 FORMAT ( //, ' HISTOGRAM OF THE MAXIMUM VELOCITY AMPLITUDE' 

1. ' DIFFERENCE DURING ACCELERATION. ', //, ' CELL SIZE', 

27X, '# OF', /, 2X, ' CCM/SEC 3 ', 6X, 'BURSTS', /) 

DO 100 I = IMIN-1, IMAX+1 
X = I - 1 
Y = X + 1. 

INUM = HIST ( I ) 

DO SO J = 1, 40 

IF( J. GT. INUM) GO TO 75 

LINE(J) = AST 

GO TO 90 

75 LINE(J) = BLNK 

00 CONTINUE 

100 WRITEC 1,101) X, Y, INUM, LINE 

101 FORMAT ( 3X, F3. 0, F3. 0, 110, 2X, ' ! ', 40A1 ) 

RETURN 

END 
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